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SECTION  I 


INTRODUCTION 


Recent  advances  in  sensor  technology  have  led  to  increased  interest  in 
the  detection  and  classification  of  various  objects  through  use  of  their 
emitted  energy.  Two  forms  of  emitted  energy  which  appear  to  be  very  promising 
for  classifying  a  wide  variety  of  objects  are  the  passively  generated  seismic 
and  acoustic  signals. 

Many  attempts  have  been  made  to  detect  and  classify  objects  by  processing 
the  signals  detected  by  geophones  or  microphones.  These  attempts  have  had 
limited  success  because  of  a  lack  of  understanding  of  the  signal  source  and 

of  the  properties  of  the  propagation  medium--in  these  cases,  the  ground  and 
atomsphere. 

Many  studies  of  pure  seismic  waves,  such  as  those  generated  by  under¬ 
ground  explosives  or  earthquakes,  have  been  conducted  and  good  seismic  models 
have  been  developed  (White,  196S) .  In  the  near  field,  however,  both  seismic 
and  acoustic  energy  can  couple  into  the  ground  and  be  detected  by  the  eeo- 
phone.  6 

There  is  an  interaction  between  the  acoustic  signal  and  the  seismic 
medium  throughout  the  length  of  the  propagation  channel.  The  signal  sources 
of  interest  generate  both  seismic  energy,  which  is  coupled  directly  into  the 
ground  and  propagated  through  the  ground  to  the  geophone,  and  acoustic  energy 
which  is  coupled  directly  into  the  ground  and  also  transmitted  through  the 
atmosphere.  Throughout  the  propagation  channel  there  is  an  interaction  be¬ 
tween  the  acoustic  signal  and  the  seismic  media  along  the  interface  between 
the  ground  and  the  atmosphere.  Part  of  the  acoustic  energy  is  absorbed,  and 
part  is  reflected.  This  interface  interaction  produces  a  nonlinear  effect, 
which  is  evidenced  by  the  appearance  of  harmonics  in  the  output  when  a  clear 
tone  is  used  as  a  stimulus. 

In  this  research,  the  Wiener  theory  of  nonlinear  system  identification 
was  used  to  characterize  the  propagation  channel  associated  with  near  field 
seismic/acoustic  propagation  and  coupling.  This  theory  requires  the  meas¬ 
urement  of  the  system  response  to  an  input  of  zero  mean  white  Gaussian  noise. 
The  system,  S,  in  this  case  is  the  ground  and  atmosphere  between  the  signal 
source  and  a  geophnne  located  distance,  r,  from  the  source.  The  input  to 
the  system  is  an  audio  system,  consisting  of  an  amplifier  and  speaker, 
driven  by  a  noise  generator. 

The  basic  mathematics  for  nonlinear  identification  theory  was  developed 
y  Norbert  Wiener  (1958).  Lee  and  Schetzen  (1965)  further  improved  the  Wiener 
theory  by  developing  a  cross-correlation  technique  which  made  it  adaptable 
to  digital  computers.  French  and  Butz  (1973)  showed  how  the  computing  time 
could  be  reduced  through  the  use  of  fast  Fourier  transforms. 
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The  potential  to  obtain  a  wide  range  of  system  response  characteristics 
from  a  single  experiment  using  nonlinear  identification  techniques  has  been 
recognized  by  researchers  in  tire  field.  Several  studies  have  centered  on 
the  identification  of  biological  systems  of  animals.  Marmarelis  (1972)  use 
the  cross-correlation  techniques  of  Lee  and  Schetzen  to  yield  a  functional 
description  of  the  horizontal,  bipolar,  and  ganglion  cells  in  the  vertebrat  * 
retina  of  the  catfish.  Marmarelis  and  McCann  (1973)  used  the  cross-correlatii  n 
techniques  to  characterize  the  visual  system  of  flies.  McCann  (19/4)  continued 
to  make  use  of  the  nonlinear  identification  theory  to  study  the  photoreceptors 
and  optokinetic  flight  response-  of  flies.  Current  research  is  underway  at  "he 
University  of  Florida  by  Brownell  and  Warren  (1976)  to  study  the  respon  e  of 
the  inner  car  of  cats  to  white  noise  stimulus. 

The  Wiener  theory  is  a  generalized  method  of  analyzing  a  nonlinear  system. 
Although  previous  use  of  the  theory  has  been  almost  exclusively  in  the  analysis 
of  biological  systems,  the  theory  is  applicable  to  most  systems  that  are  time- 
invariant  and  have  finite  memory.  Therefore,  it  covers  a  very  wide  range  o!. 
physical  systems. 

Section  IT  of  this  report,  the  background  theory  of  nonlinear  Identification 
is  briefly  reviewed.  Section  III  provides  a  description  of  the  experimental 
equipment  and  procedure.  Data  analysis  and  results  are  discussed  in  Section  IV, 
with  conclusions  and  recommendations  for  future  research  presented  in  Sections 
V  and  VI,  respectively. 
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SECTION  II 


BACKGROUND  THEORY 

Wiener  expressed  the  output  of  a  nonlinear  system  in  terms  of  an  or¬ 
thogonal  set  of  G-functionals: 


yM  =I£G"[h"'xlt)] 


(1) 


where  the  jh^J  are  the  Wiener  kernels  of  the  system  and  the  JG^J  is  a 

complete  set  of  orthogonal  functionals  when  the  input,  x(t),  is  a  white 

Gaussian  process.  Each  G  is  made  up  of  multidimensional  convolution 

n 

integrals  of  degree  n  or  less.  The  first  three  G-functionals  are: 


G 


1 


hjtTjxU-'Odx 


-Af  h2(T,T)dT  (2) 

where  A6(t)  is  the  autocorrelation  function  of  the  zero-mean  white  Gaussian 
process,  x(t),  and  6(t)  is  the  unit  impulse  function.  The  nth  degree 
G-functional  is  made  of  an  nth  order  convolution: 


hn(  Y",Tn  )  x  (t'Tl)'  *  •  x(  ^n  )dTl  1  *  *  dTn 


plus  homogeneous  functionals  of  lower  order  with  kernels  solely  dependent 
upon  h^,  derived  in  such  a  manner  that  G^  is  orthogonal  to  all  functionals 

of  degree  less  than  n. 


In  this  formulation  the  identification  problem  is  the  determination 
of  the  Wiener  kernels,  h^,  of  the  system.  The  input  and  output  of  the 
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system  must  be  clearly  defined,  and  the  input  must  be  a  zero  mean  Gaussian 
white  noise.  Note  that  if  the  input  is  not  white  but  may  be  represented  as 
the  response  of  a  linear  filter  to  white  noise,  there  is  a  procedure  dis¬ 
cussed  by  Lee  and  Schetzen  (1965)  for  whitening  the  input.  With  the  input 
and  output  known,  the  orthogonality  of  the  j G  J  allows  identification  of  the 

Wiener  kernels. 

Wiener  originally  expanded  the  kernels  in  terms  of  a  set  of  Laguerre 
polynomials,  although  any  complete  orthogonal  set  would  have  served  the  same 
purpose.  He  then  constructed  a  system  (see  Figure  1)  whose  response  to  a 
white  Gaussian  input  was  a  G-functional  with  leading  term  containing  a  kernel 
consisting  of  a  product  of  Laguerre  polynomials.  Then  by  correlating  the 
response  of  this  system  with  the  response  of  the  unknown  system  to  the  same 
input,  the  coefficients  of  the  Laguerre  expansion  could  be  determined. 


White 

Gaussian 

Process 


Figure  1.  Wiener's  Scheme  for  Measuring  Kernels 


The  identification  method  proposed  by  Wiener  allows  the  direct  measure¬ 
ment  of  the  Laguerre  coefficients  through  synthesis  by  appropriate  analog 
circuits.  Practical  considerations,  however,  require  finite  expansions 
which  result  in  unavoidable  truncation  errors.  These  errors  will  be  the 
minimum  integral  square  errors  associated  with  expansion  of  a  function  by  a 
finite  orthogonal  set.  The  Wiener  expansion  and  measurement  scheme,  there¬ 
fore,  approximates  the  kernels  in  a  minimum  integral  square  sense  over  the 
entire  domain  of  the  kernel. 

While  the  Wiener  method  is  more  applicable  to  analog  measurements,  other 
methods  have  been  devised  which  are  more  suitable  for  digital  computation 
of  the  kernels.  These  are  the  cross-correlation  method  of  Lee  and  Schetzen 
(1965)  and  the  fast  Fourier  transform  (FFT)  method  of  French  and  Butz  (1973) . 
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The  cross-correlation  method  of  Lee  and  Schetzen  allows  the  computation 
of  the  Wiener  kernels  without  expansion  in  terms  of  an  orthogonal  set.  This 
technique  may  be  used  to  compute  the  kernels  point-by-point  in  terms  of 
simple  time  delays,  o..  The  analogous  known  system  for  this  method  is  made 

up  of  parallel  time  delay  filters  (see  Figure  2)  for  which  the  impulse  response 
is  6(1-0^.  The  output  from  each  filter  is  x(t-cr). 

The  zero  order  kernel,  Gp[hg,x(t)]  =  h^,  is  the  average  value  of  y(t) 
for  the  input  x(t)  : 


h 


0 


*  y(t) 


(3) 


where  the  bar  is  used  to  denote  time  average.  Note  that  xft-o^  is  a  func¬ 
tional  of  the  first  order,  that  is 


(  t’0i)=jC  6(T‘°1  )  x  (t‘T)dT 


(4) 


Therefore ,  in  the  averaging  of  yfOxft-op,  all  averages  involving  G-func- 

tionals  of  order  2  and  higher  will  disappear  because  of  their  orthogonal 
construction.  Then 


yCOxft-Oj  )  * 


1  im  1 
T-*<»  2T 


ICjPi.xCt)] 


♦ 


dt 


1  im 

s 

T-*co 


h1 (t)A6(t-o1  )  dt 


‘  Ahi(°i) 

(5) 
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Therefore, 


h,(°l)  ■  t-Oj)  (6) 

The  second  oroer  kernel  is  computed  by  the  same  process  revealing: 

h2(v“2)  ■  7'">-h)>h)  -  v(v°2)  m 

The  second  term  in  h2  is  an  impulse  at  *  o2  which  results  from  the  av-  - 
eraging  process  of  lower  order  terms.  Since  the  lower  order  terms  are  known 

y(t)  prior  to  the  avera8ing  p™"ss  «hich  •**- 

h2(°l,02)  ■  [>,(t)-ho]x(t-0l)>(t-<,2)  (8) 

for  all  Oj,  o2. 

In  general. 


hn(°l  •  •  •  %)■  «»DS..»C«)1I»(*-01)  -XK) 

(9) 

*h?  F/FT  meth°d  °f  FIench  and  Butz  C1973)*  the  time  functions  x(t), 
y(t)  and  hn|o1  .  .  .  an j,  for  n>0,  are  expressed  in  their  Fourier  transform 

representation,  respectively,  X(f) ,  Y(f),  and  H^fj . fj  where  upper 

case  letters  are  used  to  denote  Fourier  transforms.  Also,  complex  expo- 

-j27Tf.t 

With  ira?uls!  resP°nse  e  1  ,  are  used  instead  of  the 
time  delay  filters  used  in  the  cross-correlation  method.  The  functional 
representation  of  the  output  *n(t)  of  the  known  system  then  becomes: 


r  ,  r  r  _j(fiTi  +  •  •  ♦  f  t  ) 

vt,  ^  •  •  I  .  -  <(*-'.)•  •  • 

■  x*(fj ) . . .  x^y32'^ '  ■  ■  ■' f")' 

where  the  asterisk  denotes  complex  conjugate. 


(10) 


The  zero  order  kernel  is  just  H(0)  =  Y(0)  =  y(t)  =  hn,  while  the  French 

and  Butz  xpressions  for  the  Fourier  transform  of  the  first  and  second  order 
kernels  are: 

Mj(f)  =  \  Y (f) X * f f ) 


(11) 


The  Fourier  transform  of  a  white  Gaussian  process  is  again  a  stochastic 
process;  hence,  the  kernels  calculated  by  Equation  (11)  for  a  given  sample 
input  will  represent  only  one  sample  of  the  respective  kernels.  To  obtain 
fhe  true  kernels,  an  ensemble  average  must  be  performed.  The  requirement 
for  ensemble  averaging  becomes  apparent  when  consideration  is  given  to  the 
fact  that  the  second  order  kernel  must  be  zero  for  a  linear  system.  The 
second  order  transfer  function  is  a  filter  which  passes  that  portion  of  the 
output  at  frequency  ^  ♦  f 2 )  which  is  dependent  on  the  input  at  frequencies 

f^  and  f^.  If  the  system  is  linear,  there  are  no  output  values  at  frequency 

(fj  ♦  f2)  except  those  which  are  caused  by  an  input  at  frequency  (fj  +  f2  )• 

Therefore,  the  second  order  kernel  must  be  zero  for  a  linear  system.  Without 
ensemble  averaging,  this  is  not  the  case. 


The  frequency  domain  kernels  are  therefore  defined  as 


Hj(f)  =  j  <Y(f)X*(f)> 


(12) 


and 


*  q)**(f2)> 


(13) 


The  discrete  Fourier  transform,  F(kAf),  used  in  this  computation  is 


F(kAf)  =  $  £  f(nAt)e'j2Trkn/N  (14 

n=0 

where  N  is  the  number  of  data  samples  at  sampling  interval  At  and  f(nAt) 
is  the  sampled  time  function. 


A  is  defined  as  the  power  spectral  density  (I'SD)  of  the  white  noise 
input  x(t).  The  theory  assumes  that  the  PSD  of  the  input  is  constant.  In 
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practice,  however,  the  actual  input  noise  source  exhibits  a  flat  PSD  over 
only  a  finite  range  of  frequencies  and  even  in  this  range,  the  PSD  is  not 
exactly  constant.  In  this  research  the  PSD  of  the  source  was  computed  by 
averaging  the  flat  region  of  the  sampled  PSD  from  50  to  1000  hertz  and  then 
ensemble  averaging  40  data  samples. 


When  the  kernels  and  the  input  and  output  signals  are  expressed  in  their 
respective  Fourier  transform  representation,  the  first  three  G-functionals 
from  Equation  (2)  become 


GoPVx(t)]  =  Ho 


Gjfhj.xU)] 


Hj  (f)X(f)e^2lTftdf 


dfjdf 
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SECTION  III 


EXPERIMENTAL  ARRANGEMENT 


The  experimental  arrangement  for  measuring  input  and  output  signals  is 
depicted  in  Figure  3.  The  experiment  consisted  of  playing  broadband  noise 
into  a  loudspeaker.  The  speaker  output  was  used  to  excite  the  medium.  The 
input  signal  was  generated  by  a  low  frequency  noise  generator  amplified  by 
Sansui  5000A  100-watt  amplifier.  Speaker  output  was  measured  by  a  micro- 

wereCt  rath  dlTe?tly  ™  front  of  the  sreaker.  System  output  measurements 
ere  taken  by  geophones  located  at  intervals  of  10,  20,  30,  and  50  feet  from 

the  speaker.  A  geophone  was  located  directly  beneath  the  speaker  to  measure 
seismic  energy  transferred  directly  by  the  speaker  enclosure.  Microphone 
measurements  were  also  taken  concurrently  with  each  output  geophone  for  a 
-■elated  experiment  involving  the  investigation  of  acoustic  coupling  prop¬ 
erties.  The  signals  were  recorded  on  magnetic  tape  by  an  Ampex  2200  FM 
recorder  A  Tektronix  434  Oscilloscope  and  a  Nicolet  440A  Mini  Ubiquitous 
Spectrum  Analyzer  were  used  to  monitor  the  signals  during  the  recording 
process.  b 

Ine  noise  generator  was  a  General  Radio  Company  1381  Random  Noise  Gen¬ 
erator  with  operating  modes  of  2  hertz  to  50  kilohertz,  2  hertz  to  20  kilo- 
hertz  or  2  hertz  to  2  kilohertz.  A  Tektrooics  function  generator  proved 

-  *;X^r"al  fllter  which  could  be  used  to  further  decrease  the  noise 
bandwidth. 

The  speaker  system  was  a  Pioneer  CS-99  system  which  contains  a  12-inch 
woofer  plus  a  5- inch  midrange  and  2  horn  tweeters.  The  speaker  performed 
very  poorly  w  the  lower  frequency  region.  Speaker  response  was  relatively 

anri^nn0!?  600  hertz*  but  sPeaker  power  decreased  rapidly  between  600 

and  700  hertz  due  to  crossover  effects.  The  Pioneer  speaker  was  replaced 
during  the  latter  part  of  the  test  by  a  15-inch  Altec  416-8A  speaker-.  The 
Altec  speaker  was  capable  of  responding  to  tones  as  low  as  30  hertz:  however 
a  proper  enclosure  was  not  available  for  this  speaker  at  the  time  of  the 
experiment  which  caused  poor  quality  data. 

nnm«abl!  -/r°VideS  3  c.0|"Plete  listing  of  equipment  along  with  manufacturers' 
names  and  identifying  numbers  and  the  estimated  or  rated  frequency  response 

-n  /r®iJT"ary  measuremfnts  of  tbe  system  response  to  discrete  tones  of 
30  to  2000  hertz  were  taken  to  determine  the  frequency  response  range  of  the 

Syfem  resP°nse  above  700  hertz  was  nil.  These  measurements  also 
provided  information  on  system  nonlinearity  and  for  determination  of  sampling 
rate  for  the  digital  analysis.  In  changing  from  one  tone  value  to  the  next  8 
a  frequency  sweep  was  made  to  provide  data  for  cross-correlation. 

Noise  samples  were  taken  in  the  frequency  ranges  2-2000  KHz,  2-1000  Hz, 
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impulses  created^by  stiiking'a^tfl^lSe*  with0*^'  ^  Signals  included: 
of  the  input  microphone  and  geophone ,Pen pirn,  *  hanBner  in  the  near  vicinity 
ton  van  with  enginS  «d/„  *  •«- 

signals  created  by  making  engine  on  anH  over,the  inPutsJ  and  vehicle 
speeds  by  the  inputs.  d  englne  off  passes  at  5-  and  10-MPH 

Signals  from  the  following  sensors  were  recorded  during  the  experiment: 
•The  input  from  the  noise  generator  (X). 

•Microphone  1  (Ml)  collocated  with  the  source. 

at  12*  incherbelow^ound'leveK  “i,h  th*  S°UrCe  **  Eround  level  and 

sourceMiCr°Ph°ne  2  (M2)  Space<'  »“«™«tely  at  10  and  30  feet  from  the 

source 'at  Aground  .'e^l^luo'aTu'l^chL'beloJ 
sourc,.MlCrOPh°ne  3  (M3)  S"aC<!d  alternately  at  20  and  SO  feet  from  th, 

source 'at'ground  level^nd^lso^t^l^inches^below  ground^eveT.1*1* 

Recording  arrangement  on  th,  7-ch.nnel  Amp.,  recorder  was  as  follow,: 


CHANNEL 

1 

2 

3 

4 

5 

6 
7 

Edge  Track 


SIGNAL 

Ml 

M2 

G1 

G2 

M3 

X 

G3 

Voice 


This  arrangement  maintained  the  collocate 
signals  together  on  the  sax.  tap.  h^'^  ~ 

bothD;Uh“teh”  “ere  mad,  at  each  site 

with  the  speaker  lying  flee  dom  o?thS  V**®  fr0U?d  facing  the  array  and 
and  speaker  configuration  for  each  e»P«^;J^Tlffi.1£“1*,,a 
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1  was  bordered  to  the  west  by  concrete  and  to  the  ncrth  by  an  asphalt- 

rSiSr,  wssi;  :sc^ 

“  progress  approximately  1  mile  to  the  east  which  could  be  measured  by  the 
geophone . 

Site  2  was  primarily  open  area  with  a  highway  approximately  100  feet  to 
the  ^uth  Se  »«  at  sue  2  was  primarily  sand  sparsely  covered  by  vega- 
tation.  T  e  array  was  spaced  south  to  north. 

Table  2  shows  the  site,  speaker  configuration  and  distance  from  the  source, 
and  depth  of  M  and  C3  for  each  experiment  performed,  each  experiment  eon- 
T't.d  of  the  recording  of  impulses,  clear  tones,  swept  toms,  random  noise, 

„,id  v chicle  signatures. 


TABLE  2.  DETAILED  DATA  ON  EXPERIMENTS 


i  MM  RIMENT 
;v  >IBEK 

SITE 

SPEAKER 

configuration 

G2 

DISTANCE/DEPTH 

G3 

DISTANCE/ DEPTH 

1 

1 

Upright 

10  ft/surface 

20  ft/surface 

2 

1 

Upright 

10  ft/12  in 

20  ft/12  in 

3 

1 

Facing  down 

10  ft/12  in 

20  ft/12  in 

4 

1 

Upright 

30  ft/12  in 

50  ft/12  in 

5 

1 

Facing  down 

30  ft/12  in 

50  ft/12  in 

6 

1 

Upright 

30  ft/surface 

50  ft/surface 

7* 

1 

Upright 

30  ft/8  in 

30  ft/surface 

8 

2 

Upright 

10  ft/surface 

20  ft/surface 

9 

2 

Upright 

30  ft/surface 

50  ft/surface 

10 

2 

Facing  down 

30  ft/surface 

50  ft/surface 

11 

2 

Upright 

30  ft/12  in 

50  ft/12  in 

12* 

2 

Upright 

30  ft/ 8  in 

30  ft/surface 

1 - - 

— 

J - *  — 

*Gi  located  at  30  feet  from  source,  16  inches  deep. 


The  anal  op  data  was  recorded  on  three  10-inch-diameter  reels,  designated 
„  3  Since  all  seven  channels  of  the  1/2-inch-wide  tape  were 

Tapes  -  ,  »  collection  there  was  no  room  on  the  original  tape  for  a 

Therefore^dubs  of  the  analog  data  were  made  on  1-inch-wide  tape 
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and  the  dubbed  tapes  were  time-coded.  A  strip  chart  analysis  was  nerfom,^ 
to  associate  each  signal  with  a  time  code  Each  s<k>ai  uIc  ^  rformed 

listing  of  the  signal  numbers,  the  inclusive  time  code,  and  type  signal. 


SECTION  IV 


DATA  ANALYSIS 

preliminary  analysis 

rLri°nor^'1sh:h^fdrbrry  to  dr 

™’yn£  L”c;s,‘rs,r,„r Gi  had  *•>  •“«“?* ras; 

inalyzcr  a-d  in  X-Y  plotter.  Also  thcTpsi)  ofrtheXamined  Wlt^  3  spectrum 

"oLnd  to  J“p|"8  JJ[ 

sxgnals  were  low-pass  filtered  at  1000  hertz  ?0  prevent  aliasing. 

insig^r^to^helrde'r  of^onllneaHty  1°  llTlyTel  E*  ^ 

iStS^ 

Clear  tone  responses  were  analyzed  bv  rpniavinn  ■  , 

through  a  spectrum  analyzer  and  plottingLheLesults  on  IS  x’v'JwE*8 
a  linear  scale  as  shown  in  Figures  4  through  7  At  50  hert?"  Jh  0tt?r*  Wlth 
second  harmonic  in  the  speaker  outnut  mi  8tw4  '  w  0  hertz,  there  is  a 

show  up  in  the  geophone  outputs.  This  indi^tSs^hat^heT  r CVeral  °therS 
well  as  the  seismic  medium  is  nonlinear  A  0  lnd  200 \lrrl°  *1 

harmonics  have  practically  disappeared  hut  th*  Law  h  *  ’  the  speaker 
Phones  is  still  pronounced.  At  300  Sertz  It  I  f  hai"onic  in  the  «e°- 

SL •  **?  *  • -end  harmonic^has^reappeared^in'th^speaker^^t^oo0"6 

teraction  alongLhe  interface  of^tL'cLannlLs  bSSm^^lM^walmrii  ^ 
Figure  4  indicates  that  the  nonlinearitv  nf  the  fA»0i 

—  at  so  a:  set  Alsc 
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Frequency  (.hertz) 

X  input  Signal 
Ml  Signal  at  Microphone  1 

G1  Signal  at  Geophone  1 

G2  Signal  at  feophone  2 

G3  Signal  at  Geophone  7 


Figure  5.  System  Harmonics  at  100-Hertz  Input 
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Frequency  (hertz) 
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X  Inrut  Signal 
Ml  Signal  at  Microphone  1 

III  Signal  at  Geophone  1 

G2  Signal  at  Geophone  2 

G3  Signal  at  Geophone  3 

Figure  6.  System  Harmonics  at  200-Hertz  Input 
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Frequency  (hertz) 


X 

Ml 


Input  Signal 
Signal  at  Microphone  i 

G1  Signal  at  Ceophone  1 

G2  Signal  at  Geophone  2 

Signal  at  Geophone  3 


G3 


Figure  7.  System  Harmonics  at  300-Hertz  Input 
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it  can  be  seen  from  Figure  5  that  the  third  harmonic  has  very  little  energy 
compared  to  the  second  harmonic  and  the  fundamental  at  100  hertz.  At  200 
hertz  (Tigure  6)  the  third  harmonic  in  the  geophone  signal  is  not  detectable. 
This  analysis  indicates  that  if  the  audio  system  could  be  separated  from  the 
seismic  system  there  would  be  very  little  total  energy  contribution  gained 
from  computation  of  a  third  order  kernel . 


Figure  8  is  a  plot  of  the  relative  shapes  of  the  PSD  of  the  system  signals 
taken  from  the  spectrum  analyzer.  The  speaker  output  (Ml)  is  obviously  non¬ 
white  and  therefore  does  not  qualify  as  a  proper  input  for  identifying  a  non¬ 
linear  system.  If  the  speaker  system  were  linear,  the  seismic  system  could 
still  be  identified  simply  by  using  the  input  from  the  noise  generator  and 
dividing  out  the  speaker  transfer  function,  i.e.. 


h  (f)  =  i ma 
V J  s(f) 


and 


H2(fltf2)  = 


Y(f1  ♦  f2)X*(f1)X*(f2) 


2A*S(f1)S(f2) 


(16) 


(17) 


where  Hj  and  H2  are  the  first  and  second  order  kernels  of  the  seismic  medium 

and  S(f)  is  the  speaker  transfer  function.  Since  the  speaker  system  is 
obviously  nonlinear  because  of  the  appearance  of  the  harmonics  in  the  speaker 
output,  the  above  procedure  does  not  work. 


Since  the  signal  at  G1  is  also  nonwhite  and  nonlinear,  the  same  analysis 
applies. 

Therefore,  it  must  be  concluded  that  the  transfer  characteristics  of  the 
seismic/acoustic  medium  cannot  be  determined  by  this  method  with  the  equip¬ 
ment  available.  Only  the  total  system,  including  the  audio  system  and  the 
transfer  medium,  can  be  evaluated.  The  analysis  was  continued  in  the  in¬ 
terest  of  checking  out  the  mathematical  computation  procedures  while  continuing 
a  search  for  better  equipment. 


Evaluation  of  system  response  to  impulsive  inputs  which  impart  acoustic 
as  well  as  seismic  energy  to  the  system  gives  some  indication  of  system 
memory  which,  in  turn,  indicates  the  temporal  record  length  required  to 
adequately  characterize  the  system.  Such  inputs  are  created  by  striking  a 
metal  plate  with  a  hammer.  The  acoustic  wave  which  travels  through  the 
atmosphere  is  the  first  signal  to  be  detected  by  the  geophone.  The  hammer 
blow  also  generates  a  compression  wave  which  propagates  into  the  surface  and 
a  shear  wave  which  travels  along  the  surface.  The  Rayleigh  wave  which  is 
the  vector  sum  of  the  compression  and  shear  waves  (White,  1965)  is  the  last 
disturbance  registered  by  the  geophone.  The  system  memory  is  then  defined 
as  approximately  the  time  interval  t2  -  tJf  where  t2  is  time  of  arrival  of 

the  Rayleigh  wave  and  tj  is  the  time  of  arrival  of  the  acoustic  wave. 
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Figure  8  Relative  Shape  of  Power  Spectral  Density 
of  System  Signals  and  Background  Noise 
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According  to  Marmarelis  (1972),  the  data  record  length,  R,  should  be 


R  =  MAX 


[  ■ k  •  **  *  *»] M 


(18) 


M,  the  number  of  independent  samples, 


where  f  is  the  Nyquist  sampling  rate. 

n  ^  r  &  ...  muiiucr  ui  independent  sample: 

=  brsj ts" ^ £ 

COMPUTATION  OF  THE  KERNELS 

A  first  and  second  order  kernel  was  computed  for  the  total  system  frnm 
noise  generator  to  geophone  G2  located  10  feet  from  the  speaker  * 

samples  ^The^enlels  th”  *nd  thC  ensemble  avera8e  was  tBken  over  40 

mples.  The  kernels  therefore  represent  5  seconds  of  independent  data. 

the  sr:;eps^:ru:[ng°tS:rr:LntLiinear  transfer  function  was  ~d 


L(f) 


Hf 

X( 


(19) 

where  L(f)  is  denoted  here  as  the  linear  filter  Thi*  tin.ar 

:  ss?4j.tTrU* the  abimy  °f  th;  *•  pJiis  :‘soUt. 


teiss  n  sirsrss:  smFrr 

”‘T  ?  “ITU?  oddF»r-r^" 

K  ’  G3*  G5*  and  a11  odd  order  G- functionals  contain  a 

linear  term  which  contributes  to  the  total  linear  portion  of  the  output. 

ure  if6  li*  transform  of  second  order  kernel  is  shown  in  Fig- 

scr1^vi^r'n1c^:tii:“^.r^r?:t'da^:crb;"c^idi„rt\:Ty:t'r- 

matra.s  an  the  second  order  kernel.  From  Equation  (13?  S  is  seeH  thaT 

a,  j  _ _  i  . 


H2(f2.f,). 


n.,.{ _ ...  '  -n— w*w..  is  seen  tnat 

Therefore  the  function  is  symmetrical  about  the 


Ampl i tude 
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line  fj  =  f2<  In  addition,  H 2  has  conjugate  symmetry,  i.e.,  H2(fi,f2^  = 

"f2^*  which  Allows  from  the  fact  that  is  real.  Also, 

( f i  * f 2 ^  *s  zer0  f°r  |  ^l  +  ^2  |  *  1000  because  values  for  Y(f)  are  nonzero 

only  in  the  range  |f  |  <  1000.  Using  the  above  symmetry,  H.(f  ,f  )  is  com- 

<■  X  4 

pletely  determined  by  its  values  in  Areas  1  and  2,  as  shown  in  Figure  12. 

In  order  to  evaluate  the  derived  model,  the  kernels  were  used  to  compute 
an  expected  output  for  a  given  input  signal.  The  computed  output  was  com¬ 
pared  with  the  output  measured  at  G2  for  the  same  input  signal.  Equation  (11) 
was  used  to  compute  the  output. 

At  this  point  it  is  noted  that  the  sensors  used  in  this  experiment  are 
i.ot  capable  of  measuring  a  direct  current  signal.  Therefore,  h^  *=  y^tj  is 

considered  to  be  zero  throughout  the  computation.  This  means  that  the  second 
term  of  the  second-order  kernel 

V(fi  *  f2)- 0 

and  also  the  second  term  of  the  second  order  G- functional , 

A  f  H2(f,  -f)df  -  ^ f  <Y(0)X(f)X(  -f)>df  -  0 


since  <Y (0) >  *  Hq  *  0.  The  expression  for  y(t),  the  computed  output  reduces 
to 


y(t)  =  G  (t)  ♦  G  (t) 

1  2  (20) 

where  Gj (t)  is  the  contribution  of  the  first  order  kernel  and  G2(t)  is  the 
contribution  of  the  second  order  kernel. 

The  comparisons  were  made  by  cross-correlating  the  measured  output 
signal  with  the  computed  outputs.  Figure  13  shows  the  autocorrelation  of 
the  measured  output  at  G2  compared  with  the  cross-correlations  of  the  out¬ 
put  and  that  computed  from  Hj(f)  alone,  and  then  the  output  and  that  computed 

from  both  H^(f)  and  H2(f^,f2),  Figure  14  shows  the  same  comparison  between 
the  output  computed  from  H^f),  that  computed  from  L(f),  and  the  auto¬ 
correlation  of  the  measured  signal.  A  second  comparison  of  the  outputs  can 
be  made  from  the  PSD  plots  of  Figures  15  through  18  and  Figure  19  which  com¬ 
pares  the  integrated  power  of  the  output  signals.  These  comparisons  show 


Lag  Count  (0.488  milliseconds) 


-  Autocorrelation  of  Measured  Output 

““  Cross-Correlation  of  Measured  Output  and  Output 
Computed  from 

****  Cross-Correlation  of  Measured  Output  and  Output 
Computed  from 


Figure  13.  Correlation  Coefficients  of  Computed  versus  Measured  Output 


Lag  Count  (0.488  milliseconds) 


Autocorrelation  of  Measured  Output 

Cross-Correlation  of  Measured  Output  and  Output 
Computed  from  Hj 

Cross-Correlation  of  Measured  Output  and  Output 
Computed  from  L 


Figure  14.  Comparison  of  Output  Computed  from  First  Order 
Kernel  and  Output  Computed  from  Linear  Filter 
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Figure  16.  Power  Spectral  Density  of  Output 
Computed  from  First  Order  Kernel 
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Figure  17.  Power  Spectral  Density  of  Output  Computed 
from  First  and  Second  Order  Kernels 
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Figure  18.  Power  Spectral  Density  of  Output 
Computed  from  Linear  Filter 
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Frequency  (hertz) 

Measured  Output 

Output  Computed  from  ll^f) 

“  Output  Computed  from  llj  (f)  and  ll2(f,fo) 

—  Output  Computed  from  L(f) 


Figure  19.  Integrated  Power  Spectrum  of  Measured  and  Computed  Output 


that  the  output  computed  with  the  first  and  second  order  kernel  give  a  much 
better  representation  of  the  measured  signal  than  the  output  computed  by 

be?!eenPr0CedUIT '  ,  ^  °f  the  SeCOnd  ordcr  kerncl  imPr°ves  the  correlation 

from  H  fr’alone3  rmeaSrred  °u;PUt.by  about  45  P^cent  over  that  computed 
from  Hj  Cf>  alone.  From  Figure  19,  it  is  seen  that  there  is  still  energy 

unaccounted  for  even  when  the  outputs  from  Hj (f)  and  H.,(f  f  )  are  combined. 

This  is  probably  attributable  to  the  third  and  fourth  order  nonlinearities 
which  were  observed  in  Figure  4.  "ununtarities 

COMPUTATION  PROBLEMS 

The  requirement  for  ensemble  averaging  makes  the  Fourier  transform  method 
of  computing  the  kernels  much  less  attractive  than  originally  anti ^PateS 
routing  rim,  is  increased  directly  with  the  number  of  data'sequences  u^d 
ie  ensemble.  Also,  there  is  no  good  standard  available  to  judge  the 

When *”* t h e s c S fn P , C S  rCqUired  Produce  an  acceptable  measure  of  accuracy, 
lien  these  factors  are  considered,  the  cross-correlation  method  in  time 
domain  begins  to  look  more  attractive. 

.  L^etCDr86fionmPUtinf  timeriS  exPensive*  this  is  not  the  primary  problem 
a  the  CUC  6600  computer.  The  primary  consideration  in  this  program  is 

orage  The  C DC  6600  computer  used  in  this  computation  limits  each  program 

U.  32,000  words  in  the  core  and  two  million  words  on  disk  storage  Of  course 

ti°Uld  bC  USCd’  ^  KhCn  m°re  tha"  one  taPe  is  required  to  hoZ 
ne  kernel,  the  program  must  be  interrupted  to  change  the  tapes.  This  in 

rIInn,'r,lntreaSfS  rU!"ir?  tlmC'  Routine  programs  which  exceed  10  to  15  minutes 
running  time  have  difficulty  in  being  scheduled  for  the  computer  except  at 

The  second  order  kernel  in  frequency  domain  requires  2n2  storage 

orde^kUr  W^Gre  n-1S  t5e  nUmber  °f  datU  samPles-  ^  addition,  thefirst 
order  kernel  requires  2n,  and  each  time  sequence  used  in  the  ensemble 

average  requires  n  locations.  A  half-second  of  our  data  with  40  ensemble 

averages  would  require  all  available  disk  and  core  storage  leaving  no  room 

ifrtheYem  wr  r5G  ?TOgrm’  ^is  borage  limitation  can  be  overriden 
5  P™gram  has  sufficient  priority.  A  second  possibility  is  to  reduce 
the  digitization  rate.  As  can  be  seen  from  Figure  16,  most  of  the  output 

inG^yt'S  C<?"^ln?d  tc  th?  °~  t0  300-hertz  frequency  band.  As  pointed  out 
band  tv  J1:  lnPu^n°lse  data  samples  are  available  in  the  2-  to  300-hertz 
b  JhlS  data  could  be  redigitized  at  a  600-hertz  sampling  rate  which 

would  then  better  allow  us  to  use  a  longer  data  sequence  in  each  ensemble. 
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SECTION  V 


CONCLUSIONS 


This  research  has  been  a  first  attempt  to  analyze  the  seismic/acoustic 
transfer  problem  using  the  Wiener  theory  of  nonlinear  systems.  The  research 
has  proven  to  be  very  valuable  in  that  it  has  shown  conclusively  that  the 
system  is  nonlinear,  and  linear  procedures  are  not  adequate  for  modeling  the 
system. 

The  primary  problem  area  in  this  experiment  was  the  failure  of  the  equip¬ 
ment  to  respond  in  the  frequencies  below  60  hertz.  The  targets  of  interest 
for  which  these  transfer  functions  are  required  transmit  direct  and  acous¬ 
tically  coupled  seismic  energy  of  significant  amplitude  in  the  frequency 
range  of  5  to  100  hertz  (Constantine  and  Walker,  1976). 

The  equipment  used  in  this  research  provided  no  information  below  60 
hertz.  Also  due  to  the  nonlinearity  of  the  speaker,  the  system  analyzed 
was  not  the  seismic  medium  but  the  total  system,  including  the  audio  system 
and  the  seismic  medium. 

Currently  a  search  is  underway  to  find  a  ground-vibrating  system  which 
will  respond  in  the  frequency  band  of  5  to  100  hertz.  Ideally,  this  system 
should  generate  ground  vibrations  as  well  as  acoustic  signals. 

Secondly,  a  better  means  of  measuring  the  system  input  must  be  devised. 
In  this  research  neither  the  microphone  nor  the  geophone  provided  a  complete 
measure  of  the  system  input.  The  microphone  will  not  measure  direct  ground 
vibrations  while  the  geophone,  placed  beneath  the  speaker,  is  bypassed  by 
much  of  the  acoustic  energy.  Probably  the  best  measure  of  system  input 
would  be  obtained  by  a  geophone  rigidly  attached  to  the  vibrator. 

If  the  above  equipment  problems  can  be  solved  and  a  means  of  making  a 
linear  transfer  of  energy  from  the  noise  generator  to  the  medium  devised, 
it  is  very  likely  that  the  third  and  fourth  harmonics  observed  in  this  re¬ 
search  will  disappear.  Then  the  system  can  be  adequately  characterized  by 
the  first  and  second  order  kernels  developed  here. 
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SECTION  VI 


RECOMMENDATIONS 


It  is  recommended  that: 

1.  The  effort  be  continued  to  locate  an  earth  vibrator  which  will  re¬ 
spond  in  the  frequency  range  of  5  to  100  hertz  and  make  a  linear  transfer 
of  energy  from  the  noise  generator  to  the  ground. 

2  The  computing  efforts  with  the  data  currently  on  hand  be  continued 
in  the  frequency  band  of  0-  to  300-hertz  to  refine  the  frequency  domain 
method  of  computation. 

Concurren" ly  with  frequency  domain  calculations,  a  program  should 
initiated  to  calculate  the  kernels  in  time  domain  to  provide  a  means  of 
checking  computational  efficiency,  numerical  accuracy,  and  stability  of  the 
frequency  domain  methods. 

4.  A  program  be  initiated  to  compute  the  kernels  as  a  function  of 
distance. 
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APPENDIX  A 


DATA  RECORD 


Each  signal  is  identified  by  a  5-digit  number.  The  first  digit  is  the 
number  of  the  analog  tape;  the  second  and  third  digits  represent  experiment 
number  (01  through  12) ;  the  fourth  digit  represents  signal  type  (0  for  back¬ 
ground,  1  for  sweep  time,  2  for  noise,  3  for  vehicle  data);  and  the  fifth 
digit  is  a  chronological  signal  within  each  type.  Digital  data  may  be  called 
by  signal  number.  Analog  data  may  be  located  by  its  time-coded  location  on 
the  tape. 


Tape  1 


GNAL  NO. 

TIME  CODE 

10101 

22:33:00-22:34:00 

10111 

22:40:23-22:40:33 

10112 

22:41:47-22:41:57 

10113 

22:43:50-22:44:00 

10114 

22:44:50-22:45:00 

10115 

22:46:10-22:46:20 

10116 

22:47:42-22:47:52 

10117 

22:48:57-22:49:07 

10102 

22:49:30-22:50:15 

10121 

22:50:30-22:52:30 

10122 

22:53:00-22:54:30 

10123 

22:56:00-22:57:30 

10103 

22:58:00-22:59:00 

10131 

22:59:16-22:59:36 

10132 

23:00:50-23:01:10 

10133 

23:02:13-23:02:20 

10134 

23:02: 56-23:03: 26 

10201 

23:08:00-23:08:30 

10211 

23:14:20-23:14:30 

10212 

23:15:43-23:15:53 

10213 

23:17:35-23:17:45 

10214 

23:19:07-23:19:17 

10215 

23:19:37-23:19:47 

10216 

23:21:37-23:21:47 

10202 

23:25:00-23:25:30 

10221 

23:26:00-23:27:30 

10222 

23:28:00-23:29:30 

10223 

23:29:40-23:31:00 

10224 

23:31:10-23:31:55 

10225 

23:32:40-23:33:25 

DESCRIPTION 


Background 
30-50  Hz  Sweep 
50-100  Hz  Sweep 
100-200  Hz  Sweep 
200-500  Hz  Sweep 
500-1000  Hz  Sweep 
1000-200  Hz  Sweep 
1000-100  Hz  Sweep 
Background 
0-2  KHz  Noise 
0-1  KHz  Noise 
0-500  Hz  Noise 
Background 
Vehicle  Engine  Idle 
Vehicle  Fast  Idle 
Vehicle  Power  off  Pass 
Vehicle  Power  on  Pass 
Background 
30-50  Hz  Sweep 
50-100  Hz  Sweep 
100-200  Hz  Sweep 
200- 1 00  Hz  Sweep 
100-500  Hz  Sweep 
500-1000  Hz  Sweep 
Background 
0-2  KHz  Noise 
0-500  Hz  Noise 
0-200  Hz  Noise 
0-100  Kz  Noise 
200-300  Hz.  Noise 
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SIGNAL  NO. 

TIME  CODE 

DESCRIPTION 

10301 

23:40-05-23:50:00 

Background 

10311 

23:42:22-23:42:32 

30-50  Hz  Sweep 

10312 

23:44:20-23:44:30 

50-100  Hz  Sweep 

10313 

23:47:07-23:47:17 

100-200  Hz  Sweep 

10314 

23:49:08-23:49:18 

300-500  Hz  Sweep 

10315 

23:50:25-23:50:35 

500-1000  Hz  Sweep 

10316 

23:51:53-23:52:03 

1000-700  Hz  Sweep 

10302 

23:54:00-23:54:30 

Background 

10321 

23:55:00-23:56:30 

0-2  KHz  Noise 

10322 

23:56:45-23: 57:45 

0-500  Hz  Noise 

10323 

23:58:00-23:59:30 

0-200  Hz  Noise 

10324 

00:00:00-00:01:30 

0-100  Hz  Noise 

10325 

00:02:00-00:03:30 

200-300  Hz  Noise 

10303 

00:08:00-00:08:30 

Background 

1 0331 

00:14:00-00:14:30 

Vehicle  Idle 

10332 

00:16:00-00:16:30 

Vehicle  Fast  Idle 

10333 

00: 1 7: 2 J-00: 17:30 

Engine  Off  Pass 

10334 

00:18:20-00:18:30 

Engine  On  Pass 

1 0401 

00:19:30-00:20:10 

Background 

1 0402 

00:21:30-00:22:00 

Background  (Aircraft) 

10411 

00:28:35-00:28:45 

30-50  Hz  Sweep 

10412 

00:29:30-00:29:40 

50-100  Hz  Sweep 

10413 

00:31:00-00:31:10 

100-20  Hz  Sweep 

10414 

00:34:00-00:34:10 

60-500  Hz  Sweep 

10415 

00:35:40-00:35:50 

500-1000  Hz  Sweep 

10403 

00:37:00-00:37:30 

Background 

10421 

00:43:00-00:43:30 

0-2  KHz  Noise 

10422 

00:44:30-00:45:30 

0-500  Hz  Noise 

10423 

00:46:30-00:47:15 

0-200  Hz  Noise 

10424 

00:47:45-00:48:30 

0-100  Hz  Noise 

10425 

00:48:40-00:49:30 

200-300  Hz  Noise 

10511 

00:51:15-00:51:25 

30-50  Hz  Sweep 

10512 

00:52:30-00:52:40 

50-100  Hz  Sweep 

10513 

00:55:05-00:55:15 

1000-200  Hz  Sweep 

10514 

00:56:30-00:56:40 

200-500  Hz  Sweep 

10515 

00:57:55-00:58:05 

500-1000  Hz  Sweep 

10516 

00:59:05-00:59:15 

1000-200  Hz  Sweep 

10521 

01 : 20:45-01 :21:15 

0-2  KHz  Noise 

10522 

01 : 23:20-01 : 23:40 

0-500  Hz  Noise 

10523 

01 :25:40-01 : 26: 10 

0-200  Hz  Noise 

10524 

01  -.28:00-01:28:30 

0-100  Hz  Noise 

10525 

01 : 29: 30-01: 30: 00 

200-300  Hz  Noise 

10531 

01 : 31: 30-01: 32: 00 

Vehicle  Idling 

10532 

01:34:10-01:34:20 

Vehicle  Exhaust 

10533 

01 : 34: 10-01:  34: 20 

Engine  on  Pass 

10534 

01:34:35-01:34:45 

Engine  on  Pass 

10535 

01 : 35:00-01 : 35: 1 0 

Engine  Off  Pass 

10536 

01:35:30-01:35:40 

Engine  Off  Pass 
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SIGNAL  NO. 


TIME  CODE 


DESCRIPTION 


20611 

04:27:45-04:27:55 

20612 

04:28:55-04:29:05 

20613 

04:29:10-04:29:20 

20614 

04:29:40-04:30:10 

20615 

04:31:10-04:31:30 

20616 

04:34:00-04:34:10 

20621 

04:40:30-04:42:00 

20622 

04:42:30-04:43:30 

20623 

04:44:10-04:44:40 

20624 

04:45:10-04:45:40 

20625 

04:47:00-04:48:00 

20626 

04:52:30-04:53:00 

20627 

04:53:45-04:54:15 

20628 

04:54:55-04:55:20 

20629 

04:56:30-04:57:30 

20631 

05:11:30-05:12:30 

20632 

05:13:30-05:14:00 

20633 

05:15:30-05:16:00 

20634 

05:16:40-05:17:10 

20635 

05:18:25-05:18:35 

20636 

05:19:00-05:19:10 

20637 

05:19:27-05:19:37 

20638 

05:19:53-05:20:03 

20639 

05:20:25-05:20:35 

20711 

05:22:00-05:22:10 

20712 

05:22:40-05:22:50 

20713 

05:23:30-05:23:40 

20714 

05:24:20-05:24:30 

20715 

05:31:10-05:31:20 

20721 

05:43:15-05:44:15 

20722 

05:45:00-05:46:00 

20723 

05:46:30-05:47:30 

20724 

05:48:10-05:49:10 

20725 

05:50:00-05:51:30 

20731 

05:55:13-05:55:23 

20732 

05:55:40-05:55:50 

20733 

05:56:07-05:56:17 

20734 

05:56:29-05:56:39 

20735 

05:57:05-05:57:35 

20736 

05:58:10-05:58:40 

20737 

05:59:15-05:59:45 

20738 

06:00:35-06:01:05 

20801 

06:02:00-06:03:00 

20811 

06:05:38-06:05:48 

20812 

06:09:30-06:09:40 

20813 

06:10:30-06:10:40 

20814 

06:16:20-06:16:30 

20815 

06:17:53-06:18:03 

20816 

06:18:50-05:19:10 

100-30  Hz  Sweep 
30-50  Hz  Sweep 
50-100  Hz  Sweep 
1000-100  Hz  Sweep 
200-250  Hz  Sweep 
600-700  Hz  Sweep 
0-2  KHz  Noise 
0-500  Hz  Noise 
0-200  Hz  Noise 
0-100  Hz  Noise 
200-300  Hz  Noise 
0-100  Hz  Noise 
0-200  Hz  Noise 
0-500  Hz  Noise 
0-2  KHz  Noise 
Vehicle  Idling 
Vehicle  Fast  Idle 
Vehicle  Exhaust 
Vehicle  Zooming 
Engine  Pass  5  MPH 
Engine  Out  Pass  5  MPH 
Engine  Out  Pass  10  MPH 
Engine  On  Pass  5  MPH 
Engine  On  Pass  10  MPH 
100-30  Hz  Sweep 
30-50  Hz  Sweep 
50-60  Hz  Sweep 
60-70  Hz  Sweep 
1000-100  Hz  Sweep 
0-2  KHz  Noise 
0-500  Hz  Noise 
0-200  Hz  Noise 
0-100  Hz  Noise 
200-300  Hz  Noise 
Engine  Out  Pass  5  MPH 
Engine  Out  Pass  10  MPH 
Engine  On  Pass  5  MPH 
Engine  On  Pass  10  MPH 
Engine  Idling 
Engine  Zooming 
Engine  Exhaust 
Engine  Exhaust  Zooming 
Background 
30-50  Hz  Sweep 
50-30  Hz  Sweep 
30-100  Hz  Sweep 
100-200  Hz.  Sweep 
200-300  Hz  Sweep 
1000-300  Hz  Sweep 
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SIGNAL  NO. 


TIME  CODE 


DESCRIPTION 


20817 

20818 
20819 
20821 
20822 

20823 

20824 

20825 

20831 

20832 

20833 

20834 

20835 

20836 

20837 


06:19:47-06: 19:57 
06:21 : 38-06:21 : 43 
06:22:30-06:22:40 
06:39:30-06:41:00 
06:42:00-06:42:30 
06:44:30-06:45:00 
06:46:30-06:47:00 
06:49:00-06:49:30 
06:51:45-06:52:15 
06:52:30-06:53:00 
06:53:50-06:54:10 
06:54:20- 06:54:40 
06:58:14-06:58:24 
06:59:30-06:59:46 
07:01 :00-07:02:00 


300-500  Hz  Sweep 
500-700  Hz  Sweep 
700-1000  Hz  Sweep 
0-2  KHz  Noise 
0-500  Hz  Noise 
0-200  Hz  Noise 
0-100  Hz  Noise 
200-300  Hz  Noise 
Engine  Idle 
Engine  Zooming 
Engine  Exhaust  Idle 
Engine  Exhaust  Zooming 
Engine  Out  Pass 
Engine  On  Pass 
Engine  On  Pass 


SIGNAL 


TIME  CODE 


DESCRIPTION 


30901 

30911 

3091 2 

309 1 3 

30914 

30915 

30916 

30917 

30918 

30921 

30922 

30923 

30924 

30931 

30932 

30933 

30934 

30935 
30926 
30937 

31001 

31011 

31012 

31 0 1 3 

31014 

31015 

31002 

31021 

31022 


17:46:00-17:47:00 
17:48:38-17:48:48 
17:51:50-17:52:00 
17:55:30-17:55:40 
17:59:00-17:59:10 
18:01:52-18:02:02 
18:06:20-18:06:30 
18:08:50-18:09:00 
18:10:30-18:10:40 
18:16:30-18:18:00 
18:20:45-18:21:30 
18:22:00-18:22:30 
18:28:25-18:28:55 
18:32:25-18:32:35 
18:33:40-18:33:50 
18:34:20-18:34:30 
18:35:50-18:36:00 
18:37:20-18:37:30 
18:39:14-18:39:24 
18:40:59-18:41:09 
18:43:30-18:44:30 
18:46:55-18:47:05 
18:49:40-18:49:50 
18:51:36-18:51:46 
18:53:30-18:53:40 
18:56:25  18:56:35 
19:01:10-19:02:10 
19:07:10-19:07:40 
19:09:00-19:09:30 


Background 
30-40  Hz  Sweep 
50-70  Hz  Sweep 
90-100  Hz  Sweep 
100-150  Hz  Sweep 
200-300  Hz  Sweep 
1000-400  Hz  Sweep 
400-500  Hz  Sweep 
700-1000  Hz  Sweep 
0-2  KHz  Noise 
0-500  Hz  Noise 
0-200  Hz  Noise 
200-300  Hz  Noise 
Engine  Idle 
Engine  Zooming 
Engine  Fast  Idle 
Engine  Exhaust  Idle 
Engine  Exhaust  Zooming 
Engine  On  Pass  5  MPH 
Engine  On  Pass  5  MPH 
Background 
30-50  Hz  Sweep 
50-70  Hz  Sweep 
80-50  Hz  Sweep 
90-100  Hz  Sweep 
200-300  Hz  Sweep 
Background 
0-2  KHz  Noise 
0-500  Hz  Noise 


SIGNAL  NO. 


TIME  CODE 


DESCRIPTION 


31023 

31024 

31111 

31112 

31113 

31114 

31115 

31116 

31117 

31121 

31122 

31123 

31124 

31131 

31132 

31133 

31134 

31135 

31136 

31137 

31138 

31221 

31222 

31223 

31224 

31225 

31231 

31232 

31233 

31234 


19:11:00-19:11:30 
19:26:35-19:27:05 
20:22:00-20:22:10 
20:22:50-20:23:00 
20:23:35-20:23:45 
20:26:45-20:26:45 
20:28:48-20:28:58 
20:30:47-20:30:57 
20:31:40-20:31:50 
20:35:35-20:36:15 
20:37:50-20:38:20 
20:38:55-20:39:25 
20:40:00-20:40:30 
20:42:30-20:42:40 
20:42:55-20:43:05 
20:44:00-20:44:10 
20:44:25-20:44:35 
20:45:22-20:45:32 
20:46:07-20:46:17 
20:47:03-20:47:13 
20:47: 53-20:48:03 
20:48:40-20:49:40 
20:50:30-20:51:30 
20:52:20-20:52:50 
20:53:10-20:53:40 
20:54:05-20:54:35 
20:55:00-20:55:10 
20: 55:35-20:55:45 
20:56:02-20:56:12 
20:56:40-20:56:50 


0-200  Hz  Noise 
200-300  Hz  Noise 
30-40  Hz  Sweep 
40-50  Hz  Sweep 
50-70  Hz  Sweep 
90-100  Hz  Sweep 
400-200  Hz  Sweep 
400-500  Hz  Sweep 
700-800  Hz  Sweep 
0-2  KHz  Noise 
0-500  Hz  Noise 
0-200  Hz  Noise 
0-100  Hz  Noise 
Engine  Idle 
Engine  Zooming 
Engine  Exhaust  Idle 
Engine  Exhaust  Zooming 
Engine  Out  Pass  5  MPH 
Engine  Out  Pass  10  MPH 
Engine  On  Pass  5  MPH 
Engine  On  Pass  10  MPH 
0-2  KHz  Noise 
0-1  KHz  Noise 
0-500  Hz  Noise 
0-200  Hz  Noise 
0-100  Hz  Noise 
Engine  Idle 
Engine  Zooming 
Engine  Out  Pass 
Engine  On  Pass 
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INITIAL  DISTRIBUTION 


Hq  USAF/RDQRM 
Hq  USAi'/SAMI 
Hq  USAr/XOXFM 
AFIS/INTA 
AFSC/DLCA 
AFSC/IGFG 

afsc/sdza 

afal/dho 

AFWAL/TECH  LIB/FL2802 

ASD/ENFEA 

FTD/PDXA-2 

AFML/LTM 

AFWL/NSC 

AFWL/NSE 

AFWL/SUL 

AUL  ( AU/LSE- 70-239 ) 

DDC 

OGDEN  ALC/MMWM 

Plcatlnny  Arsenal/SARPA-TS-S#59 
Army  Mat  Sys  Anal  Agcy/DRXSY-J 
Army  Mat  Sys  Anal  Agcy/DRXSY-A 
tones  Scl  Info  Ctr/Ch,  Doc  Sec 
USAE  Waterways  Exp  Stn/WESFE 
NRL/Code  2627 

NavAIr  Sys  Comd/Tech  Llb/AIR-954 
NavSurWpnCen/Tech  Lib  &  Info 
Svs  Br 

NavOrdStn/Tech  Lib 
NavA1rTestCen/CT-176  TID/Tech 
Pubs 

USNWC  (Code  533) 

Sandla  Lab/Tech  Lib  Dlv  3141 
The  Rand  Corp/Llbr-D 
TACTEC,  Battelle  Col  Lab 
TAWC/TRADOCLO 
AFATL: 

DL 

DLOSL 

DLY 

DLJ 

DLJC 

DLJF 

DLJK 

DLJM 


2 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

1 

1 

1 

2 

1 

1 

1 

2 

1 

2 

1 

1 

1 

1 

1 

1 

9 

1 

1 

1 

1 

1 

1 
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(The  reverse 


